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ABSTRACT 

One of the well-known problems of producing instruments for Extremely Large Telescopes is that their 
size (and hence cost) scales rapidly with telescope aperture. To try to break this relation alternative new 
technologies have been proposed, such as the use of the Integrated Photonic Spectrograph (IPS). Due to 
their diffraction limited nature the IPS is claimed to defeat the harsh scaling law applying to conventional 
instruments. The problem with astronomical applications is that unlike conventional photonics, they are not 
usually fed by diffraction limited sources. This means in order to retain throughput and spatial information 
the IPS will require multiple Arrayed Waveguide Gratings (AWGs) and a photonic lantern. We investigate 
the implications of these extra components on the size of the instrument. We also investigate the potential 
size advantage of using an IPS as opposed to conventional monolithic optics. To do this, we have constructed 
toy models of IPS and conventional image sliced spectrographs to calculate the relative instrument sizes and 
their requirements in terms of numbers of detector pixels. Using these models we can quantify the relative 
size/cost advantage for different types of instrument, by varying different parameters e.g. multiplex gain and 
spectral resolution. This is accompanied by an assessment of the uncertainties in these predictions, which 
may prove crucial for the planning of future instrumentation for highly-multiplexed spectroscopy. 

1. INTRODUCTION 

Spectroscopy is the most useful technique in modern Astronomy, addressing problems ranging from funda- 
mental cosmology to exoplanets. In its simplest form a dispersive spectrograph contains four components. 
A slit to isolate the area to be dispersed, a collimator, a grating or prism to disperse the light and a camera 
and detector to record the intensity at each wavelength. In order to retain throughput the slitwidth of a 
spectrograph is usually matched to the FWHM of the seeing. If the system is not diffraction limited then 
the collimated beam must increase in size in order to maintain the same resolution.^ This means that the 
disperser must also grow in size. This leads physical problems such as stresses and flexure in the materials, 
along with the difficulties inherent in building these large monolithic structures. In order to use the same 
design principles as existing instruments more exotic materials and construction processes are needed, which 
drives the costs of building the instruments up. This leads to the cost of instruments scaling with the square 
of the telescope aperture or more.^ 

Many ways to break this relationship have been considering including image slicing^ and adopting technolo- 
gies developed by the photonics industry.^ The solution considered in this paper makes use of the Arrayed 
Waveguide Gratings (AWGs) as the dispersive element. These form part of the Integrated Photonic Spec- 
trograph (IPS) concept.^ 

The IPS takes light from an input fibre which is usually matched to the seeing limit (as with conventional 
fibre fed instruments) and so supports many modes. As the input to the AWG must be single moded the 
light has to be split into individual single moded fibres by a photonic lantern.^ These are then fed into the 
AWGs which disperse the light into individual spectra.-^ 
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While the principles have been demonstrated so far these have only used single or few modes from a 
single input, so have suffered from loss of throughput and coverage of field. 

In this paper we shall look at whether photonic spectroscopy can compete with conventional instruments 
in terms of size. We will construct a toy model of conventional instruments and photonic ones in order to 
evaluate the relative volumes produced by these instruments. 



2. INSTRUMENT SIZE SCALE RELATIONSHIPS 

It has been stated in Astrophotonics papers that using AWGs breaks the telescope diameter, spectral re- 
solving power relation.^ The relation. 



xDt ~ xDt ^ ' 

holds for a long slit spectrograph using a diffraction grating as the dispersive element. Here R is the 
resolution of the instrument, m is the diffraction order, A is the wavelength, W is the length of intersection 
between the grating and collimated beam, x is the angular slitwidth, Dt is the diameter of the telescope, 
7 is the blaze angle and Dcoi the diameter of the collimated beam. It shows that for a given resolution, 
angular slitwidth, and blaze angle, that as the diameter of telescope increases the diameter of the collimated 
beam must increase, leading to a bigger instrument. Unlike a conventional spectrograph the AWG must be 
diffraction limited (A = x^t) due to its single mode nature, so the spectral resolving power can be shown 
to be of the form 



(2) 

where N is the number of waveguides, m is the diffraction order and C errors from manufacturing.^ It 
must be noted that this is for a single AWG. When telescope is not at the diffraction limit multiple AWGs 
(and hence modes) would be needed to conserve throughput. 

The number of spatial modes, in a stepped index fibre can be approximated as 



T/2 

M = ^^. (3) 



The associated V parameter being given by 



Vfibre = — . (4) 

Where s is the diameter of the fibre, and B the numerical aperture at which the fibre is operated.'^ This 
must be less than the maximum numerical aperture of the fibre. This can be modified so that the variables 
so they are same as Eq. ^ by taking 



e 



1 

= Xfr = xFtDt 



(5) 
(6) 
(7) 



with Ft being the focal ratio, /t the focal length of the telescope and s the diameter of the fibre (equivalent 
to the slitwidth in Eq. ([T])), yielding 



M^i^-^l (8) 

Therefore it can be seen that for each samphng element the number of modes increases as the square of 
the telescope diameter in a similar way to the number of slices at the diffraction limit (x^t/1-22A)^. This 
means unless we will need to utilise more AWGs in order to maintain throughput, therefore the size of the 
IPS will still scale with telescope diameter. 



3. MODELLING 

We model both the standard instrument and its AWG equivalent. For our standard instrument we use a 
double pass echelle model similar to that developed in Ref. [3l We choose a variable focal length for the 
camera (and hence telescope) to keep the model simple and compact. For our AWG model we choose a 
reflective design similar to those described in Refs. ITOlandlTTI 

The major difference between the models is the way they are scaled. All other parameters are based upon 
existing instruments for both models. As such extra components such as the optics, detector and housing 
are accounted for in the echelle model simply by using scaling factors. The instrument itself will have been 
optimised for size and structural integrity so it will represent the optimal instrument design. 
By comparison there are no full scale photonic spectrographs currently in existence. The current generation 
are prototypes, only sampling a very small area of the sky and not making use of enough AWGs to allow 
high throughput. Because of this we use a bottoms up approach when scaling our AWG model. We use 
an existing AWG obtained from Gemfire Livingstone and use its known size to scale our model. This does 
however mean the AWG model does not include the extra equipment that would be required for a full 
instrument, such as detectors and the housing of the instrument. We also assume all the AWGs in the final 
device can be stacked on top of each other, which in practise may not be possible. 

3.1 The echelle Model 

The echelle model is adapted from that shown in Ref. [121 Here we will be calculating the model using 
the initial input parameters of Ac, AXfsr^ p, where AXfsr is the free spectral range (FSR). As 

stated in Ref. 12 modelling a spectrograph for monochromatic light is simple, but for a large wavelength 
range we need to add in extra parameters accounting for dispersion. For our model here we use ±^/2 to 
represent the extra beamspread due to the wavelength dispersion and ±^/2 to represent the dispersion due 
to conservation of etendue. Our toy model is based on an echelle spectrograph with double pass construction 
(see Fig. [1]). For an unsliced spectrograph its dimensions are estimated as 



Lx = {Dcoi + ^ {ciech + Dcoi tan 7) + 

^ech) ^ech 

(9) 

Ly = {Dcoi + ^ {aech + Dcoi tan 7) + aech) bech (10) 
Lz = {DcoiFc^ Dcoi tsin J ^2aech) bech' (11) 

Where Dcoi is the diameter of the collimated beam, 7 the blaze angle and Fc the focal ratio of the collimator 
and camera. We include aech and bech as an oversizing parameter in order to allow the model to be fitted 
to existing instruments, with aech as the additive term and bech the multiplicative. We consider Littrow 
configuration at blaze, so we can re-arrange Eq. ([T]) to give the diameter of the collimated beam 

Dcoi = (12) 
'^''^ 2 tan 7 ^ ^ 

The blaze angle in Littrow configuration is sin 7 = , and the FSR of a system is AXfsr = — • Thus 



T = tan7 = tan (sm~^ ( — • 



(13) 



To work out the spatial beamspread we use conservation of etendue, yielding 



^ = NyX^. (14) 

J^col 

We now work out <l>, using the FSR and Ac- In order to calculate the total angle occupied on the detector 
by the spectrum we use the grating equation sin(^^) + sin(^o) = mpX. Setting Oi = j gives 



^ = 0{Xmax)-0{Xrmn)- (15) 

Where Xmax and Xmin are the maximum and minimum wavelengths that the system operates at (= Ac ± 
Aa/2). In order to find the minimum focal length for our system we need to make sure we can resolve each 
point in our spectra. We know that x^cDt = nodp^ where no is the oversampling of the detector (we use 
2.5 allowing for a non-gaussian shape in the point spread function (PSF)) and dp is the size of each pixel 
(our default is 13.5/im). From this we can calculate the minimum focal ratio of the camera Fc. 

With all this in place we can rewrite the equations for the scale lengths of the instrument as 



LAx) = (^^ + ^{a,ch + ^^^+aechyech (16) 

Lyix) = I +^laee/t + j+aechjbech (17) 

Lzix) = (^^^ + l) + 2aee/.) 6ec.. (18) 



Finally to correctly sample the output we require one sampling point in the spatial direction (as we have 
fixed our slit the the FWHM already) and no samples in our spectral direction, yielding 



AT AT AT ^>^FSRno AXpsuRnQ . . 

Npixeis noNyNx — . (19) 

ZAA ^min 



3.2 The AWG Model 



AWGs are usually purpose built and their design is hugely varied, especially when it comes to arrangement 
allowing the path difference between waveguides. To keep the toy model simple we have chosen a refiective 
AWG with Rowlands Circle arrangement for the Free Propagation Region (FPR).^^'^ We have removed the 
bend at the refiective end of the waveguide array for simplicity (see Fig. [2]). Because of this it looks almost 
identical to a double pass echelle spectrograph model described above. 

Using the definitions in Fig. [2] we arrive at the following equations for the size of the model AWG 



Lx = {max{D, E) + aawg) bawg (20) 

Ly = {Cawg + '^)bawg (^1) 
Lz = {aawg^A^{N^g-l)AL)bawg^ (22) 



Here N^g is the number of waveguides, A is the x length of the FPR, AL the length difference between 
adjacent waveguides to achieve the required order for a given X^ D is the length containing the waveguides 



(analogous to the illuminated length of the grating in the echelle model), E the x length of the detecting 
surface, w is the waveguide diameter. The oversizing parameters are aawg^ ^awg and Cawg- 

First we calculate the appropriate dispersion order we use the re-arranged FSR equation for an AWG 

(23) 



^^FSR 

Here m is the diffraction order and Uc is the refractive index of the waveguides. This can be combined 
with Eq. (j2j) and D = N/p yielding 

D=^. (24) 

mp 

We can make use of the equation for free spatial range in an arrayed waveguide grating {Xfsr = 
XminLpp/ns), where is the refractive index of the FPR and Lf the length of the free space propagation 
region. Combining with geometrical arguments gives 

E = Lf sin (^^^^y (25) 
In order to calculate AL we make use of the equation for calculating the central wavelength of the AWG 

AL ^— (26) 

Tie 

Ac the central operating wavelength. A can be calculated using geometry and is 

A = LfCos{e)=LfCos(^-^y (27) 

Where N is calculated using Eq. (j2j). In order to calculate Lf we make use of the fact the imaging 
requires the number of detector pixels to be able to adequately sample the PSF at the resolution required 
(equivalent to sampling of the echelle model in Ref. 3). To do this we take the dispersion relation 



SX 



(28) 



^5x J LfTTlp 

This can be combined with Eq. (|23|) to remove m. We take the equation for resolution SX = Xmin/R and 
Sx = nodp. , We use Xmin in order to get the maximum Lf This allows Eq. (|28|) to become 

UsUodpRAXpsR 
> -T2 • (29) 

r min 

Finally we can calculate the number of pixels we need for the required resolution 

Npixels = = ' (30) 



4. RESULTS 

We set the input parameters in Table [1] to represent GNIRS pausing those hsted in Ref. '3^ as a guide. We 
have selected GNIRS as it is representative of an existing facility-class instrument. We then vary resolution, 
field of view (FOV) of the instrument, telescope diameter and seeing. 

The results are presented in graphical form in Figs. [3]to[T0l The immediate conclusion is that in all cases 
the AWG instrument will be smaller. This is slightly misleading though as we have not included detector 
volumes or instrument housing, which would increase the scale lengths, leading to a much closer comparison. 
In all cases the number of detector pixels is much higher for the AWG model, though if the input is close to 
the diffraction limit this number decreases, making the devices more comparable. 

Figure [3] shows varying the resolving power will produce similar scale length differences between the AWG 
and echelle models. This means there is no significant advantage for using the AWG instrument for low or 
high resolution spectroscopy. In Fig. |4]we see that using instruments with higher resolving power also results 
in a similar difference between the AWG and echelle model, suggesting that resolution will not be a major 
factor in deciding the use of the instrument. 

Figure [5] suggests the size of the AWG instrument will have the greatest advantage for smaller fields of 
view. The reason for this is that for small fields of view the angular dependence of the echelle makes the 
overall instrument larger than the AWG. When the slit length is increased the increase in the angle for the 
echelle becomes smaller than the extra size of the extra AWGs. This is illustrated in Fig. [TTl Figure [6] shows 
once again the number of detector pixels on the AWG instrument becomes larger as the FOV increases, but 
this is at the same rate as the echelle instrument. 

Figure [71 suggests increasing the size of the telescope will not make any significant difference to the size 
difference of the AWG and echelle models. This means the same size issues will be present for AWG instru- 
ments as conventional ones for ELTs using the model setup. We also see from Fig. [8] that as we increase the 
telescope diameter the number of pixels required increases with the number of propagating modes (see Eq. 

©)• 

Figure [9] suggest that the relative sizes of the instruments is not highly dependent on seeing. This should 
be qualified though. If the input is at the diffraction limit there will be no need for a photonic lantern, or 
potentially a smaller one, leading to a smaller and less complicated instrument, if the input is far from the 
diffraction limit it will require a large photonic lantern. Figure [10] shows the number of detector pixels is 
comparable for diffraction limited inputs, but increases hugely as the seeing becomes worse, due to the extra 
modes (see [8]). 



5. CONCLUSION 

We have used simplified models for both a conventional echelle model and an AWG model to simulate the 
size and number of detector pixels of GNIRS on Gemini. We have then varied resolution, telescope diameter, 
sample size and seeing. Our results are presented in terms of the scale length (cube root of the volume) of 
the instrument and the number of detector pixels the instrument will require. 

We find that the AWG model will be slightly smaller in size for all cases becoming much smaller when 
the FOV of the instrument is small. This is slightly misleading though as we have not included other 
components in the model (such as detectors and housing). This means the actual instrument will most likely 
be the same size if not larger than a conventional non-photonic counterpart. We also find this problem will 
be compounded for larger telescopes unless the input is very close to the diffraction limit. We also find 



that larger a larger FOV makes the AWG instrument less competitive in terms of size, meaning that AWG 
instruments may be suited to instruments with smaller fields of view, such as amateur instruments or for 
high resolution spectroscopy. 

Our results indicate the number of detector pixels will be substantially more for photonic instruments unless 
the input is very close to the diffraction limit. This suggests that either near diffraction limited Adaptive 
Optics or a space based mission would be required for optimal performance from photonic spectrographs. 

In conclusion, our findings suggest that photonic spectroscopy would be most useful on small instruments 
with diffraction limited seeing. We find that the AWG instrument will be on the same magnitude in terms 
of size as a conventional spectrograph, but will have more detector pixels. 

This does not provide a comprehensive model as we have only studied photonic spectrographs with a single 
input. Recent Astrophotonics papers'^ ^ suggest that multiple inputs could be used with varying configura- 
tions. This would potentially reduce the size of the instrument, but would not reduce the number of detector 
pixels required unless all the spectra could be combined onto a single detector pixel. We intend to investigate 
this further with similar models in the near future. 
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Figure 1. The echelle model. The left part of the image shows the spectra dispersion of the echelle model, and the 
right part shows the spatial extent of the device. This is identical to the one described in Ref. [121 




Figure 2. The AWG model, Left is the x-y view of the AWG, the top right shows the side view of the AWG, with the 
waveguide in the center. Bottom right is an enlargement of the FPR. 
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Table 1. Table of input parameters for GNIRS calibration, taken from Ref. "T2"and calculated using the model. These 
will be used as the default parameters for the models. 
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Figure 3. The resulting scale lengths due to varying 
the resolution of the GNIRS models with the different 
scaling parameters. This shows the AWG instrument 
will be smaller over all resolutions. The dashed line 
represents GNIRS. 




Figure 4. The resulting pixel numbers due to varying 
the resolution of the GNIRS models with the dif- 
ferent scaling parameters. This shows the AWG in- 
strument requires orders of magnitude more detector 
pixels. The dashed line represents GNIRS. 
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Figure 5. The resulting scale lengths due to varying 
the slit length of the model instruments. From this 
we can see that increasing the slit length will have a 
far greater effect on the AWG model than the con- 
ventional instruments. The dashed line represents 
GNIRS. 
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Figure 6. The resulting pixel numbers due to vary- 
ing the slit length of our model instruments. From 
this we can see that increasing the length will have 
a far greater effect on the AWG model than the con- 
ventional instruments. The dashed line represents 
GNIRS. 




Figure 7. The resulting scale lengths due to varying 
the diameter of the telescope feeding our model in- 
strument. The difference in instrument sizes is not 
greatly influenced by telescope diameter. The dashed 
line represents GNIRS. 



Figure 8. The resulting pixel numbers due to varying 
the diameter of the telescope feeding our model in- 
strument. The different in number of pixels increases 
greatly when the telescope diameter increases. The 
dashed line represents GNIRS. 
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Figure 9. The resulting scale lengths due to varying 
the seeing of the telescope. The advantage of using an 
AWG instrument is greatest when the input is close 
to the diffraction limit. The dashed line represents 
GNIRS. 
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Figure 10. The resulting pixel numbers due to vary- 
ing the seeing of the telescope. It is easily seen that 
unless the input is near the diffraction limit the num- 
ber of pixels in the system will be much greater in 
the AWG device. The dashed line represents GNIRS. 
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Figure 11. An example of changing the sht length of the different model instruments, using a slit length of Uy yields 
a larger echelle instrument, though increasing the slit length by a factor of 5 makes the AWG instrument larger. 



